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ABSTRACT 

We present a Gaussianity analysis of the WMAP 5-year Cosmic Microwave Back- 
ground (CMB) temperature anisotropy data maps. We use several third order estima- 
tors based on the spherical Mexican hat wavelet. We impose constraints on the local 
non-linear coupling parameter fni using well motivated non- Gaussian simulations. We 
analyse the WMAP maps at resolution of 6.9 arcmin for the Q, V, and W frequency 
bands. We use the KQ75 mask recommended by the WMAP team which masks out 
28% of the sky. The wavelet coefficients are evaluated at 10 different scales from 6.9 to 
150 arcmin. With these coefficients we compute the third order estimators which are 
used to perform a x^ analysis. The x^ statistic is used to test the Gaussianity of the 
WMAP data as well as to constrain the fni parameter. Our results indicate that the 
WMAP data are compatible with the Gaussian simulations, and the fni parameter is 
constrained to —8 < fni < +111 at 95% CL for the combined V+W map. This value 
has been corrected for the presence of undetected point sources, which add a positive 
contributio n of A fni = 3 + 5 in th e V+W map. Our results are very similar to those 
obtained by iKomatsu et al.l (|2008l ) using the bispectrum. 

Key words: methods: data analysis - cosmic microwave background 



^ 



1 INTRODUCTION 

A valuable source of information of the early universe is 
the cosmic microwave background (CMB) radiation. The 
temperature fluctuations of this radiation can be used 
to put tight constraints on the cosmological parameters. 
The Big Bang theory and the inflationary models are 
our best theories to describe the universe. I n part i cular , 
the standard, single field, slow roll inf l ation (|Guthl llQSlI : 
lAlbrecht fc SteinhardtllTosJ : iLinddllQSl Il983h is one of the 
most accepted models because of the accuracy of its pre- 
dictions. In particular, it forecasts that the primordial den- 
sity fluctuations are compatible with a nearly Gaussian ran- 
dom field. However, different primordial processes (as for ex- 
ample topological defects, etc.) can introduce no n- Gaussian 
featu res at certain levels that may be detected (Cruz e t al.l 
120071 1. Also, several non-standard models of inflation predict 
a detectable level of non- Gau ssianity in the prim ordial grav- 
itational potential (see, e.g. iBartolo et al.l 12004 and refs. 



therein). Other non-Gaussian deviations can be explained 
by the presence of foreground contamination or systematic 
errors. In any case, the search of non-Gaussian deviations in 
the CMB has become a question of considerable interest, as 
it can be used to discriminate different possible scenarios of 
the early universe and also to study the secondary sources 
of non- Gaussianity. 

Many tests of Gaussianity have been performed on dif- 
ferent CMB data sets. For a recent review of this topic see for 
example Mart inez- Gonzalez (2008). Some of the tests have 
found that the data are compatible with Gaussianity but 
there are some important detections of non- Gaussian devia- 
tions (see for exam ple Eriksen et al 



2004, 



2005; Cop i et al 
2004 . "20 061: 1 Vielva et al II2004I: IC ruz et al. 20051: IWiaux et al. 



2006 : .Vielva et al.ll2007l : iMonteserin et al- 2008l). all of them 



using the Wilkinson Microwave Anisotropy Prob^ (WMAP) 
data. 

High precision experiments as the WMAP are sen- 
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sitive to deviations due to second-order effects in per- 
turbation theory, usually parametrised through the lo- 
cal non-Unear coupling parameter fni. Several methods 
have been used to constrain this parameter using the 
data of different experiments. We can mention the an- 
gular bispectrum on WMAP data ( Koniatsu et aU l2003l: 
Creminelh et all l2QQ6l: ISpergel et"ai1 boOTl : iKomatsu etaJl 
2008) and COBE JKomatsu et all l20Q2l): the Minkowski 
functional on WMAP (Ko matsu et a l. 2003; Spergel e t al. 
2007 : Gott et al. 2007 : Hikage et al.1 2008 : Komatsu et al 



20081). BOOMERa n G (IPe Troia et al.ll2007l ) and Archeops 
(jCurto et a l* '2007", '2008*); differe nt kind of wavelet s on 
WMAP (^ukheriee & Wang 200l ICabella et aDl2005l ) and 
COBE (Cavon et al.l 120031 ) among others. The previous 
works suggest that fni is compatible with zero at least 

at 95% confidence l e vel. However oth er recent works 

(| Jeong fc SmootI l2007l : lYadav fc Wandeltl l2008l l suggest a 
positive detection of fni at a confidence level greater than 
95% using WMAP data. 

In this paper we perform a wavelet-based analysis of the 
5-year WMAP data in order to constrain the fni parameter. 
We use the high resolution WMAP data maps (6.9 arcmin) 
and realistic non- Gaussian s imulations performed follo wing 
the algorithms developed bv lLiguori et al.1 (|2003l . l2007l ) . 

Our article is organised as follows. Section [2] presents 
the statistical method, the estimators that we use to test 
Gaussianity and constrain fni , as well as the data maps and 
simulations to be analysed. In Section [3] we summarise the 
results of this work. The conclusions are drawn in Section |4l 



2 METHODOLOGY 

2.1 The SMHW 

For this analysis we use a wavelet-based technique. 
The considered wavelet is the Spherical Mexican Hat 
Wavelet (SMHW) as defined in iMartmez- Gonzalez et al.l 
((2002). The spherical wavelets have been used in some 
analyses to test the Gaussianity of different data sets. 
We can mention the analysis of the COBE-DMR data 
(iBarreiro et all l2000l: ICavon et al.1 l200lL [2003) and WMAP 



data ([Vielva et al 



Cruz et al. 



2004 l2007l: iMukheriee fc Wand 



Cabella et al.ll2005l: ICavon et al.ll2005l: iMcE wen et al.ll200^ 
20051 . I2OO7I : IWiaux et al.1 [2OO8) among oth- 



200 



ers. The SMHW can be obtained from the Mexican hat 
wavelet in the plane E? through a stereographic projection 
(JAntoine fc VanderghevnstI 119981 1. Given a function /(n) 
evaluated on the sphere at a direction n and a continuous 
wavelet family on that space ^(n;b,i?), we define the con- 
tinuous wavelet transform as 



/ 



w{h-R) = / dn/(n)*(n;b,7?) 



(1) 



where b is the position on the sky at which the wavelet co- 
efficient is evaluated and R is the scale of the wavelet. In the 
case of the SMHW we have that th e wavelet only depends on 
the p olar angle 6 and the scale R ((Martinez- Gonzalez et al.l 
I2OO2I I. In particular we have 



^5 



-R) 



27tN{R) 



-(1)1 h-m 



-yV2R^ 



(2) 



N{R) = R 



and 



y^2tan(^-j 



, R' R^ 



1/2 



(3) 



(4) 



We can compute the wavelet coefficient map given by Eq. \T\ 
at several different scales Rj in order to enhance the non- 
Gaussian features dominant at a given scale. In particular, 
we win consider 9 scales: Ri = 6.9', R2 = 10.3', R3 = 13.7', 
R4 = 25', R5 = 32', Re = 50', R7 = 75', Rs = 100', and 
Rg = 150'. In addition we will consider the unconvolved 
map, which will be represented hereafter as scale Rq. Larger 
scales are less sensitive to the local fni model. 

2.2 The estimators 

In previous works, constraints on fni with wavelets have 
been obtained using only the sk ewness of the wavelet co- 
efficients at different scales (e.g. IMukheriee fc W ang 2004 
■Cabella et al.,,2005, ). In the present work, other third order 
moments involving different scales are also considereqj. As 
it will be shown in the results section, the combination of 
these estimators is as efficient as the bispectrum. The esti- 
mators that we use in this analysis are based on third order 
combinations of the wavelet coefficient maps Wi{Rj) evalu- 
ated in sets of three contiguous scales. For each scale Rj and 
the next two scales Rj-\-i and Rj-\-2 we can define 

N ■ —1 

q2{Rj,Rj+i) = 2^ — i- 

i—O ■^ 

N ■ —1 

qsiRj^Rj+i) = ^T^ V 
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1 
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-1 . 



crjCrj+i(Tj+2 



WijWi,j+2 
J- 3+2 



qQ{Rj,Rj+2) 



N^ 



1 v;^ ^^i,j^^ 

,3+2 ^ cr?cr 

N ■ - 

- E 



N ■ —1 

■'■^ptx -*- 2 

W^,jWiJ^2 



3,3 + 2 



3^j+2 



(5) 



where 



where Nj^j^j^ is the number of available pixels after com- 
bining the scales i?j^ , Rj^ , and Rjr^ , Npix is the total num- 
ber of pixels, Wij = Wi{Rj), and aj is the dispersion of Wij. 
Each map Wij is masked out with an appropriate mask at 
the scale Rj and its mean value outside the mask is removed. 
These estimators have a Gaussian-like distribution when are 
computed for a set of Gaussian simulations. Thus, we can 



^ Notice that inter-scale wavelet estimators have been previ- 
ously used in the context of blind Ga ussianity analyses, see e.g. 
iPando et all (| 19981 ): IMukheriee et al.' f2000) (for analyses involv- 
ing two scales) and [Cav on et al. (2001) (involving three scales). 



© 0000 RAS, MNRAS 000, 000-000 



WMAP 5-year constraints on fni with wavelets 3 



use effectively a x^ statistics to test Gaussianity and to con- 
strain fni . Considering all the estimators evaluated in all the 
scales we can construct a vector 

V = {qi{Ro),q2{Ro,Ri),q3{Ro, i?i), , q^Ro, i?i, i?2), •.•) (6) 

with a dimension of Uv — rise + 2 (rise — 1) + 3 (rise — 2) for 
Use ^ 2, where Usc is the number of considered scales. This 
vector is used to compute the x^ estimator 



riv-l 



X^ = ^ {vk - {vk))Cf^i\vi - {Vl)) 



(7) 



where () is the expected value for the Gaussian case and Cki 
is the covariance matrix Cki = (vkVi) — {vk){vi). The Gaus- 
sianity analysis consists on computing the x^ statistic for 
the data and compare it to the distribution of this quantity 
obtained from Gaussian simulations. 

The second part of the analysis consists on setting con- 
straints on the fni through a x^ test. In this case 



X^ifnl) = ^(Vk- {Vk)Ui)Cj^i\fnl){vi - {Vl)f^,) 



(8) 



where ()/^^ is the expected value for a model with fni and 
Cki = {vkVi)f^i - {vk)f^i Mur Fo^ low values of fni {fni ^ 
1500) we can use the following approximation Cki{fni) — 
Ckiifni = 0) = Cki- The best-fit fni for the data is obtained 
by minimization of X^ifni)- Error bars for this parameter 
at different confidence levels are computed using Gaussian 
simulations. 



2.3 Data and simulations 

For our analysis, we use the 5-year WMAP foreground re- 
duced data, which are available in the LAMBDA web sitqj. 
We combine the maps of different radiometers, using the in- 
ve rse of the noise variance as an optimal weight, as described 
in Be nnett et al.l (|2003|). In particular we use five combined 
maps: Q+V+W, V+W, Q, V, and W. The Q+V+W map 
has the eight maps corresponding to the two radiometers of 
the Q band (41 GHz), two radiometers of the V band (61 
GHz) and four radiometers of the W band (94 GHz). Sim- 
ilar combinations are constructed for the V+W, Q, V and 
W maps. The pixel resolution of these maps is 6 .9 arcmin, 
corresponding to a HEALPix (JGorski et al.ll2005r ) Nside pa- 
rameter of 512. The mask that we use is the KQ75 which 
discards 28% of the sky. 

We should take into account the effect of the mask in the 
wavelet coefficient maps. Pixels near the border of the mask 
(the Galactic cut and other features) are affected by the zero 
value of the mask (|Vielva et al.ll2004l l . Therefore we need one 
extended mask for each scale that removes these affected 
pixels from the analysis. We use the method described in 
iMcEwen et al._ (200 5) to construct our extended masks. We 
compute the wavelet coefficients at each scale of the KQ75 
mask without the holes corresponding to the point sources, 
and consider only the pixels with low values (that is, the 
less affected pixels). Then we multiply by KQ75 to mask 
out the point sources. In particular, our threshold is 0.001, 



http://lambda.gsfc.nasa.gov 



i.e., all the pixels which have a wavelet coefficient for the 
KQ75 mask larger than 0.001 in absolute value are masked 
out. In addition we also test the effect of the extension of the 
mask by applying a less restrictive threshold of 0.01. In Fig. 
[1] we present the masks that we use for the 10 considered 
scales. 

We also need Gaussian and non- Gaussian simulations 
for this analysis. The Gaussian simulations are performed 
as follows. For a given power spectrum d (we use the best 
fit power spectrum for WMAP provided by LAMBDA), we 
generate a set of Gaussian aim- From these multipoles we 
produce a map for each different radiometer by convolving 
with the corresponding beam transfer function. We also in- 
clude the pixel properties by convolving with the pixel trans- 
fer function. We add a Gaussian noise realisation to each ra- 
diometer simulation and then we combine them in the same 
way as the data maps. Following the analysis of the WMAP 
team, we assume that the instrumental noise is well approx- 
imated by Gaussian white noise at each pixel. This noise is 
characterised by a dispersion that depends on the pixel po- 
sition and the corresponding radiometer. Although the data 
also contain small residuals of 1// noise (IJarosik et al.ll2003l , 
2007) their contribution here is expected to be negligible. 

The non- Gaussian simulations are produced following a 
model which introduce s a quadratic term in the primordial 
gravitational potential (ISalopek fc Bond! 199C)l:lGangui et al.l 



Il994l : IVerde et al.ll2000l : iKomatsu k Spergelll200lf h 

$(X) = $^(X) + /n^{$i(x) - ($i(x))} (9) 

where $l(x) is a linear random field which is Gaussian dis- 
tributed and has zero mean. This kind of non- Gaussianity 
is generated i n various non-stand ard infiationary scenar- 
ios (see, e.g. iBartolo et al.l l2004r ). The simulations with 
fni are gene r ated f ollowi ng the algorithms described in 
iLiguori et al.l (|2003l , l2007h . In particular we have a set of 

(G) 

300 Gaussian simulations in the a^rn space, a)^ and their 
corresponding non-Gaussian part a^^ . A simulation with 
a given value of fni is constructed as 



.(G) 



+ fnia 



(NG) 



(10) 



To produce each non-Gaussian WMAP simulation with fni 
we transform the multipoles defined in Eq. [10] into each dif- 
ferencing assembly map using the corresponding beam and 
pixel transfer functions and we add a Gaussian noise reali- 
sation. Then, we combine these maps to form the Q+V+W, 
V+W, Q, V and W combined maps. 



3 RESULTS 

In this section we present the Gaussianity analysis of the 
WMAP data using the Q+V+W, V+W, Q, V, and W com- 
bined data maps. First we consider the estimators defined 
in Eqs. \S\ and the x^ statistic given in Eq. [71 and compare 
the values obtained for the WMAP data with the distri- 
bution constructed with Gaussian simulations. Second we 
constrain the local non-Unear coupling parameter fni using 
realistic non- Gaussian simulations. Finally we estimate the 
contribution of the undetected point sources to the best fit 
value of fni. 
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Rq = 50' 



R^ = 75' 




R„ = 100' 



Rp = 150' 




Figure 1. Masks considered for the analysis at each wavelet scale. The case Rq is the KQ75 mask. The others correspond to the 
restrictive extended masks (black) and the less restrictive extended masks (white+black). The grid corresponds to a size of 20 degrees. 



3.1 Analysis of WMAP data 

We estimate the wavelet coefficient maps at the 10 consid- 
ered scales defined in Subsect 12121 In particular, for Usc = 10 
we have 52 estimators as the ones defined in Eq. [5] In Fig. 
[2] we present these statistics for the WMAP data compared 
with the values obtained for 1000 Gaussian simulations. We 
analyse the Q+V+W, V+W, Q, V, and W combined maps 
using the restrictive extended masks, and we also analyse the 



V+W map with the less restrictive extended masks (marked 
in the tables and figures with {V + W)*). The qi estimator 
(defined in Eq. \5\ corresponding to the skewness) has a sim- 
ilar shape th an the one obtained for the 1-year a nd 3-year 
data maps (IVielva et al.1 l2004l : iMukheriee &: WansL 200i : 
ICruz et aLll2007f ). In all the cases the skewness obtained for 
the data is compatible with the skewness of Gaussian simu- 
lations. We can see in Fig.[2]that the values of the estimators 
with k ^ 36 are systematically below the mean. However, for 
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Figure 2. The 52 statistics Vk, given in Eq. [G] evaluated at 10 different scales for the WMAP combined maps. From left to right and top 
to bottom, we present the values corresponding to the Q+V+W, V+W, V+W, Q, V, W maps using the restrictive extended masks, and 
the V+W map using the less restrictive extended masks. The diamonds correspond to the data. We also plot the acceptance intervals 
for the 68% (inner in red), 95% (middle in green), and 99% (outer in magenta) significance levels given by 1000 Gaussian simulations of 
signal and noise. 



a correct interpretation of this trend it is important to take 
into account that those estimators are strongly correlated at 
large scales. In particular, considering the normalised corre- 
lation matrix in the V+W case, dj /{aiaj), we have that 
its elements are dj /{aiGj) ^ 0.7 on average for vector in- 
dices greater or equal than 36, whereas dj /{aiGj) ^ ^A for 
vector indices lower than 36. 

We have also performed a x^ analysis in order to check 
if the data are compatible with Gaussian simulations. This is 
plotted in Fig. [3] We estimate the mean and the covariance 
used in Eq. [71 with another set of 1000 Gaussian simulations. 
We can see that for all the cases, the data are compatible 
with the Gaussian simulations. The statistical properties of 
the histograms are presented in Table [1] where we give the 
X^ value obtained for the data, the degrees of freedom of the 
X^ statistic (DOE), the mean and the dispersion obtained 
for Gaussian simulations and the cumulative probability for 
the data. 

Finally we impose constraints on the fni parameter us- 
ing the considered maps. We calculate the expected values 



of all the estimators for different fni caseqj (see left panel of 
Fig.Ufor the combined V+W map). Then, we perform a x^ 
analysis using Eq.[8]in order to find the best fit value for fni 
for each data map. We also analyse Gaussian simulations 
in order to obtain the frequent ist error bars. In the central 
and right panels of Fig. |4]we present a plot of X^ifni) ver- 
sus fni for the combined V+W data and a histogram of the 
best fit fni for 1000 Gaussian simulations of the combined 
V+W map. Table [2] lists the fni values which best fits the 
five considered maps, and the main properties of the distri- 



"^ We use 300 simulations to calculate the mean values of the 
estimators for different fni cases and 1000 Gaussian simulations 
to calculate the covariance matrix. We analysed the convergence 
of the mean values by computing the best-fit fni value of the data 
using two independent sets of 150 simulations. The difference in 
the obtained value is lower than 1 and therefore even with ~100 
simulations the convergence is achieved. Analogously, for the co- 
variance matrix we have checked that convergence is achieved 
with ~1000 Gaussian simulations. 
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Figure 3. The distribution of the x^ statistics obtained from Gaussian simulations of the Q+V+W, V+W, Q, V, W maps using the 
restrictive extended masks, and the V+W map using the less restrictive extended masks. The vertical lines correspond to the values 
obtained from the data. 



but ions of best-fit fni obtained from Gaussian simulations. 
From this table we have at 95% CL that —8 < fni < 118 for 
the combined V+W map, -42 < fni < +93 for Q+V+W, 
-63 < fni < +102 for Q, -46 < fni < +109 for V, and 
—38 < fni < +114 for W (using the restrictive extended 
masks defined in Subsect. l2.3p. These values are compatible 
with the ones obtained bv iKomatsu et alj (|2008l ) using the 
bispectrum. Notice that fni increases as we go from Q to 
W maps. The fni va lue obtained for the V +W map is the 
same as the one that iKomatsu et al.l (|2008) obtains for this 
map using imax = 700 and the KQ75 mask. The Q+V+W 
map has a smaller fni value. The best-fit fni increases with 
frequency for the Q, V, and W maps, and the a{fni) of the 
Q map is larger than the one o f the V and W maps. This 
also agrees with IKomatsu et al.l (|2008 ) . 

The less restrictive extended masks give tighter results 
for the constraints on fni compared with the restrictive ones. 
In particular, from Tabled we have that —5 < fni < +114 
at 95% CL. In this case the dispersion of fni is smaller 
since the available area is larger for all the scales. The ad- 
ditional pixels are not significantly affected by the mask for 
this analysis and the wavelet is still efficient. These con- 



strain ts are very similar to those obtained with the bispec- 
trum ([Komatsu et al.ll2008l l, which are —5 < fni < +115 at 
95% CL. 

Notice however that other works, e.g. lYadav fc Wandeltl 



(2008), have found evidences of non-zero value for fni at 2.8cr 
CL using the KSW bispectrum. We do not find deviations 
with respect to the zero value using the 5- year data at 95 
% CL and the same result is presented in Komatsu et alJ 
(|2008l V We may wonder if that deviation is introduced by 
the 3-year WMAP mask (KpO), uncertainties in the char- 
acterisation of the beams, or some systematics present in 
the 3-year data. To check this possibility we analysed the 
V+W 3- year WMAP data (updated after the 5-year release) 
using our estimators, the KpO mask and its corresponding 
extended masks. The results are — 17 < fni < +108 at 95% 
CL, which are compatible with the 5-year results. Therefore 
our method is not able t o detect the deviations reported in 
lYada v fc Wa ndeltl (|200 8^1 using the 3- year data. This may be 
explained by the fact that the wavelet method probes dif- 
ferent combinations of scales than the bispectrum and also 
responds differently to possible systematics present in the 
data. 
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Figure 4. From left to right, the normahsed mean values of the estimators for 300 simulations of the V+W map with different f^i 
contributions, the x^ifnl) statistics for the WMAP combined V+W data map, and the histogram of the best fit fni values for a set of 
1000 Gaussian simulations of the V+W map. We have used the less restrictive extended masks. Similar results are obtained with the 
restrictive extended masks. 

Table 1. x^ constructed from the Sstatistics for Q+V+W, V+W, Q, V, and W data maps using Eq. [7] We also present the mean and 
the dispersion of the x^ corresponding to 1000 Gaussian simulations, and the cumulative probability for the x^ of the data. 
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^data 


DOF 


ix') 


a 


P{X 


^^xLJ 


Q+V+W 


54.65 


52 


55.49 


11.78 




0.63 


V+W 


43.23 


52 


55.22 


11.73 




0.20 


Q 


62.84 


52 


55.31 


12.08 




0.86 


V 


49.81 


52 


55.24 


11.96 




0.44 


w 


41.74 


52 


54.92 


11.34 




0.16 


(V+W)* 


55.31 


52 


55.47 


12.47 




0.65 



The V+W map is analysed using the less restrictive extended masks. 



We have also tested other third order estimators that 
involve the derivatives (in particular the squared modulus 
of the gradient and the Laplacian) of the wavelet coefficient 
maps. These estimators are less sensitive to the fni param- 
eter and its error bars are significantly wider than the ones 
obtained with the estimators of Eq. \5\ 



3.2 Contribution of undetected point sources 

It has been shown that the Mexican hat wavelet is 
very useful in de tecting point sources (Cavo n et al.1 12000 : 
IVielva et al.ll200lh since its contribution is enhanced in the 
wavelet coefficient maps at certain scales. Therefore it is 
important to study the contribution of the undetected point 
sources on the fni estimation. The skewness and kurtosis 
of wavelet coefficient ma p s due to point sources have been 
studied in lArgiieso et al.1 (120061) . We estimate the contribu- 



tion of the point sources to our estimators by performing 
Monte Carlo simulations. We use a straightforward model 
to simulate the radio sources following the one described in 
iKomatsu et al.l (|2008[ ). In this model it is assumed that all 
the sources have the same intensity flux, {Fsrc =0.5 Jy), 
and that they are randomly distributed in the sky follow- 
ing a Poisson distribution. The source contribution to the 
temperature at any pixel is given by 



Tsrc{n) 



ih(x/2)' 



24.8 MJy/K Qp 



Att/N^ 



ptx 1 



(11) 



e is a Poisson 



where x — i//(56.8 GHz), Qpix - 
random variable of mean (e) = Qp 
sity of sources. The power spectrum for this type of source 



and Usrc is the den- 



distribution can be easily calculated (|Tegmark fc Efstathioul 
Il996l ). 






,2sinh^(x/2) 



67.55 MJy 



(12) 



Using risrc = 85 sr~^ and a pixel resolution of Nside =512 
we have a power spectrum of Cf^ — 8.68 xlO~^ jiK^ sr 
for the Q band. This model roughly reproduces the mea- 
sured values of the power spe ctrum and the bispectrum 
(|Nolta et aDl2QQ8l : lKomatsu et al. 2008 ). 

For this analysis we have performed 1000 point source 
simulations for the V+W map. We add these point source 
simulations to the V+W simulations of CMB with noise. For 
each one we estimate its best fit fni and compare it with the 
values obtained for the same simulation but without point 
sources. The difference gives us an estimate of the contam- 
ination due to the point sources. They add a contribution 
of A/ni = 11 ± 4 to the CMB considering the restrictive 
extended masks. For the less restrictive extended masks the 
point source contribution is A/^z = 3 =b 4. The difference 
with the previous result is explained because the zero val- 
ues of the KQ75 mask affect the efficiency of the wavelet 
to detect point sources. This can be seen in the left panel 
of Fig. [5] where we plot the mean values for the third order 
estimators for the V+W simulations including point sources 
for both kinds of extended masks. The mean value of sev- 
eral estimators, specially the skewness for small scales, is 
lower for the less restrictive extended masks. This reduces 
the Afni value. Considering the contribution of the point 
sources, our estimate of fni for this case is +22 < fni < +83 
at 68% CL and -8 < fni < +111 at 95% CL. For compari- 
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Table 2. Best fit f^i values obtained from Q+V+W, V+W, Q, V and W combined maps. We also present the mean, dispersion and 
some percentiles of the distribution of the best fit fni values obtained from Gaussian simulations. 



Map 


best fni 


{ful) 


<yUni) 


^0.160 


^0.840 


^0.025 


^0.975 


Q+V+W 


26. 


-2. 


34. 


-34. 


30. 


-68. 


67. 


V+W 


58. 


-1. 


33. 


-35. 


34. 


-66. 


60. 


Q 


16. 


-1. 


42. 


-43. 


41. 


-79. 


86. 


V 


30. 


0. 


39. 


-39. 


38. 


-76. 


79. 


w 


40. 


-4. 


39. 


-44. 


35. 


-78. 


74. 


(V+W)* 


56. 


0. 


30. 


-31. 


30. 


-61. 


58. 



The V+W map is analysed using the less restrictive extended masks. 



son, the bias intr oduced by point sour ces for the bispectrum 
is A/ni = 5 =b 2 (|Komatsuet al.ll2008i) . 

We also estimate the contribution of undetected point 
sources using the more realist ic source number c ounts 
dN/dS, derived from the work of lde Zotti et alJ (|200 5Y The 
dependence of dN/dS with the frequency is very small in 
the range between 61 GHz and 94 GHz. Thus, for simplic- 
ity, we assume the same dN/dS for both V and W maps, 
evaluated at a frequency of 71 GHz. We select a range of in- 
tensities between Smin — 1 mJy and Smax — 1 Jy- Then we 
generate point sourc e simulations followin g the distribution 
dN/dS derived from Ide Zotti et all (|2005l ). The simulations 
are transformed from intensity flux to temperature as in the 
previous case (following Eq. Ill|) . For this analysis we have 
also performed 1000 point source simulations for the V + VK 
map and have added them to Gaussian CMB plus noise 
simulations. We analyse the resultant simulations as in the 
previous case to obtain the contribution to fni due to the 
point sources. They add a contribution of Afni = 3 =b 5 to 
the CMB considering the less restrictive extended masks and 
A/ni = 17 it 5 considering the restrictive extended masks. 
As in the previous case, the differences in Afni using the 
restrictive and less restrictive extended masks are explained 
because the zero values of the KQ75 mask affect the effi- 
ciency of the wavelet. This is plotted in the right panel of 
Fig. El 

Therefore, considering the realistic point source model 
and the less restrictive masks our estimate of fni remains 
unchanged (-8 < fni < +111 at 95% CL). 



4 CONCLUSIONS 

We have analysed the 5-year WMAP data using Gaussian 
and realistic non- Gaussian simulations through a wavelet- 
based test. We have considered different combined maps for 
the analysis: Q+V+W, V+W, Q, V, and W. We have con- 
sidered two kinds of extended masks to analyse the V+W 
map (defined in Subsect. l2.3|) . The third order moments (Eq. 
O of the wavelet coefficients of these maps are compatible 
with Gaussian simulations (see Fig [2]). 

We have performed a x^ analysis and found that the 
data are indeed compatible with Gaussian simulations (see 
Fig-EJ- We performed another x^ analysis to constrain the 
non-linear coupUng parameter fni by using non-Gaussian 
simulations with fni- The best-fit fni values of the five 
analysed maps are co mpatible with the ones obtained by 
iKomatsu et al.l (|2008r ) using the bispectrum, showing sim- 



ilar confidence intervals and also a similar trend with the 
frequency. 

Finally we have estimated the contribution to fni from 
unresolved point sources for the V+W map using a sim- 
ple model that has point sources with con stant intensity 
Fsrc = 0.5 Jy) and a realistic model given bv lde Zotti et al.l 
20051 ) . We have found that they add a positive contribution 
of Afni = 11 ±4 for the simple model and Afni = 17di 5 for 
the realist ic model. The se values are larger than the one ob- 
tained bv lKomatsu et al.l (|2008[ ) for the bispectrum, which 
can be explained by the enhancement of the point sources 
produced by wavelets. Using the less restrictive extended 
masks, the point source distribution add a positive contribu- 
tion of Afni = 3±4 for the simple model and Afni = 3=b5 for 
the reahstic model. The smaller values are explained because 
the less restrictive extended masks add some pixels that af- 
fect the efficiency of the wavelet to detect point sources. 
Taking into account the point source correction, our best 
estimate of fni is — 8 < fni < +111 at 95% CL. It is im- 
portant to emphasise the agreement found between the two 
estimators (bispectrum and wavelets) , since they are formed 
by very different combinations of the data and are affected 
by systematic effects (like the mask, noise, beam response 
and foreground residuals) in very different ways. 
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Figure 5. Normalised mean value of the third order estimators for 1000 Gaussian simulations plus point source simulations of the V+W 
combined map using the restrictive mask (threshold of 0.001) and the less restrictive mask (threshold of 0.010). In the left panel we 
use the model with point sources of constant flux, Fsrc=0.5 Jy, and in the right panel we use the model with a dN/dS derived from 
Ide Zotti et al.l (|2005h . 
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